A GAUGE-INDEPENDENT MECHANISM FOR CONFINEMENT AND MASS GAP: 

PART II — G = SU{2) AND D = 3 



J. WADE CHERRINGTONi 



^Department of Mathematics, University of Western Ontario, London, Ontario, Canada 

Abstract. We apply to the case of gauge group G = SU{2) in three dimensions a recently proposed gauge- 
independent mechanism for confinement that is based on a particular form of the dual spin foam framework 
for lattice gauge theory. Explicit formulae for interaction factors and their asymptotics are introduced and 
their behavior in different sectors of the theory are identified and analyzed. We arrive at several elementary 
properties of the dual theory that represent one scenario by which confinement may be realized at weak 
coupling. We conclude with an outlook for further development of this approach. 



1. Introduction 

The present work provides a detailed application of a new framework for investigating the confinement 
of static charges and mass gap in Yang-Mills theory to the case of gauge group G — SU{2) and space- 
time dimension £) = 3. This framework was presented in [7] and is based on an exact duality between 
conventional Yang-Mills in the lattice regularization and a lattice spin foam model. The case of G = SU{2) 
in three dimensions is particularly attractive as a first testing ground for the general case, as explicit vertex 
amplitudes are known and possess a particularly tractable form At the same time, this case is rich 

enough to exhibit in a non-trivial fashion the various features introduced in [7]- A previous application of 
spin foams to the confinement problem for this case has also been made based on analogy with the dual 
photon and monopole picture of confinement in U{1) theories 8 ; the present work also uses the spin foam 
transformation but proceeds in a different direction. 

Broadly speaking, the results we obtain here can be summarized as follows. Using an exact duality in 
terms of spin foams we show explicitly from first principles how an area law behavior for the expectation 
values of Wilson loop observables may be obtained in the weak coupling limit. While short of a rigorous 
proof, we can clearly identify properties that would be required to provide such a proof — they are given 
as conjectures. Thus, from the perspective of finding a rigorous proof of confinement, what we have done is 
present a new framework, within which can be formed a reduction of the confinement problem into several 
"smaller" conjectures that are elementary propositions regarding the spin foam model. In doing so, we 
hope to invite attention to this new perspective that stands apart from many past approaches in offering an 
understanding of confinement and mass gap that is manifestly gauge- invariant. 

The structure of the paper is as follows. In Section 2, we recall the spin foam expansion of the expecta- 
tion value of a Wilson loop observable. Applying the general framework of [7^, the worldsheet-background 
decomposition and interaction factors are found explicitly. In Section 3, we analyze the interaction factors 
in more detail, and in particular introduce an asymptotic formula for the vertex interaction Iv{f) that holds 
in the large spin limit. Applying asymptotic estimates for the vertex interaction, we find that for the class 
of self-avoiding worldsheets, one can show that /\/ — > 1 in the weak coupling limit — a property we shall 
refer to as "vertex decoupling" , a special case of the vertex stability property introduced in [7] ■ In Section 
4, we consider different sectors of the full theory according to the framework described in [7] . In the case 
of self-avoiding worldsheets interacting with asymptotic background, we show how area-law behavior is par- 
ticularly transparent due to vertex decoupling. The remaining sectors of the theory are also discussed and 
specific properties are identified as necessary conditions for confinement. In Section 5 we give formal state- 
ments of these required properties in a form suitable for pursuing via analytic methods or through numerical 
method such as the recently developed non-abelian dual simulation algorithms J5] . Section 6 concludes with 
an assessment of what has been shown and discusses ways in which the results may be strengthened. 



2. The lattice spin foam form of the Wilson loop observable 

We start by introducing the passage from the conventional to spin foam form of LGT, and review the 
worldsheet-background formulation of the Wilson loop observable that was introduced in [7]. 

2.1. The spin foam dual. As reviewed in [7J and references therein, the spin foam transformation of lattice 
Yang-Mills theory can be summarized as follows: 
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where gp is the holonomy consisting of a product of group valued variables and their inverses along the edges 
of p. In the present case of G = SU{2) and £> = 3, we have 
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The first line of ([T]) is the conventional lattice regularization of a gauge theory using a discretized action S{gp). 
The second line represents the character expansion of the action that has often been used in strong coupling 
expansions for lattice gauge theory [10] . The final line is based on a lattice spin foam duality, which can be 
understood as a refinement of the character expansion method; Consistent with the formalism introduced 
in [7j, for a given spin foam / the notation f{p), /(e), and f{v) denotes the set of irrep and intertwiner 
variables local to a plaquette, edge, and vertex respectively. A general description of this transformation can 
be found in |20| while details of the transformation for our present case is discussed in ^-5,9 . 

The configuration space of the spin foam model is Fy , the set of admissible spin foams — all colorings 
of the plaquettes and edges by unitary irreps of 5/7(2) satisfying certain local constraints. Admissibility 
constraints are necessary but not sufhcient conditions for the edge and vertex amplitudes (and hence overall 
amplitude) to be non- vanishing; a discussion can be found for example in [Sj. 

To provide a concrete model, we assume an intertwiner basis has been chosen, in which case the above 
diagrams become spin network evaluations that can be efficiently computed on the computer, as described 
in [5]. Here and throughout we use the character expansion coefficients Cj^ from the heat kernel action [15j . 

A commonly used order parameter for confinement (and the one we shall use in what follows) is the 
expectation value of a Wilson loop observable. This expectation value is defined as follows: 
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where Oesr ^5 (Se) denotes a product of representation matrices (g) in the ^th (fundamental) irreducible 
representation of 

As shown in [SlU] , the presence of field insertions associated to the Wilson loop leads to an expansion in 
spin foams 
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where the admissibility conditions of spin and intertwiner labels are changed along the Wilson line and 
modified edge and vertex amplitudes and Ae and Ay are to be used. The barred vertex amplitude has 
the form given in Figure [T] for edges and vertices in F, with the vertex amplitude depending on whether 





Figure 1. Charged vertex amplitudes: examples of perpendicular and parallel threading. 

the Wilson loop enters and leaves in the same direction or makes a perpendicular turn at the vertex. For 
the barred edge amplitude along F, an insertion of a line of half- integer charge modifies the vacuum edge 
amplitude given in ([3]). For edges and vertices not in F, the barred amplitudes are equal to the original 
unbarred amplitudes. In summary, we have rewritten the conventional form of the Wilson loop expectation 
value ([5]) as a ratio of the partition function of two different systems: 
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Here we note that Zy is built of up of only vacuum configurations Fy, and Zr is built up of configurations 
Fr that satisfy modified admissibility conditions (note that Fy and Fr are disjoint sets). Since Zy is 
independent of F, the dependence of the expectation value on F is controlled by Zr (a given lattice and fixed 
/3)- ^ 

As shown in [7j , the spin foam transformation offers us a view of the Wilson loop expectation value that 
shares many similarities with the strong coupling expansion but retains its validity in the continuum limit. 
We describe this next. 

2.2. Definitions and decomposition of F-spin foams. Let Fr denote the set of all spin foams admissible 
in the presence of the given Wilson loop F [T- admissible spin foams), and let Fy denote the set of spin foams 
admissible in the case of no external charge; we shall refer to these as vacuum admissible spin foams. Finally, 
let Sr denote the set of spin foams which can be viewed as fundamentally charged (possibly self-intersecting) 
worldsheets ending on F; we shall refer to these as T-attached worldsheets. 
We next introduce an important fact that was shown in Theorem 3.1 of [7]: 



Wilson loops in higher representations can also be analyzed by the methods discussed here but won't be detailed in the 
present work. 

3 



Fact 2.1. Worldsheet-background decomposition. Every element of Fr can be expressed as the sum of 
a vacuum spin foam and a F-attached worldsheet. That is, there exists a map p from Fr into Fy x Sr-' 

(8) p:Fr^FvxSr 
such that the inverse map Fr satisfies V(/, s) G p{Ft) 

(9) P''[if,s)] = f + s, 

where the set p{Fy) is the image of the set Fr under p and the + operation is the addition of half-integer 
labels of a configuration from Sr at the corresponding plaquettes and intertwiner labels at the corresponding 
edges of the lattice k. 

We shall refer to the p in the above as a "decomposition" . Defining a decomposition is equivalent to 
specifying a pair of functions we shall call f{-y) and 5(7) that pick out the Fy and Sr components of 7 G Fr, 
respectively. 

We should emphasize here that as discussed in [7] , there is no unique way of choosing this decomposition. 
A given decomposition p = (/(7), 3(7)) constitutes one (of many possible) maps with the desired properties. 
As discussed in [7|, one way of dealing with this degeneracy is to carry out unrestricted sums over Fy x Sr, 
but weight every configuration by a factor accounting for the degeneracy. This is the strategy we will apply 
below; others may be possible. 

2.3. The Worldsheet-background factored amplitude. By making use of the above result, we can 
rewrite Zp as a double sum over vacuum admissible spin foams Fy (which in the present context we shall 
refer to as "background" spin foams) and the set of F-connected worldsheets Sr as follows: 
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Here we have expressed the expectation value of a Wilson loop F in terms of a double sum: a sum over 
worldsheets whose boundary is compatible with F and a sum over background spin foams (from the vacuum 
ensemble). Because the sum over Sr is unrestricted, the same foam from Fr may come up multiple times: 
we define this degeneracy to be D{f,s). We thus divide by this degeneracy factor to avoid over-counting; 
see Appendix A of [7] for further discussion of degeneracies. 

Observe that the presence of the worldsheet modifies the background amplitude by a factor that is localized 
to the worldsheet. This factor we shall refer to as the interaction factor, it represents the re- weighting of 
the background amplitude due to the presence of the worldsheet. We can now write the expectation value 
of the Wilson loop as an expectation value on the dual ensemble as follows: 



(12) (OrJ = 
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where we define the E observable (a function on the dual spin foam ensemble Fy) as the sum all worldsheets 
bounding F of the interaction factor weighted by the inverse degeneracy factor: 
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We next examine the interaction factor (and in particular its weak coupling limit) in more detail. 
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3. Analysis of Interaction Factors 



In this section we take a closer look at the interaction factors coming from the plaquette, edge, and vertex 
amplitudes. In particular, we shall present the large spin behavior of these quantities that play a decisive 
role at weak coupling. 

3.1. Plaquette interaction. In what follows we will use a "twice-spins" convention of using integers rather 
than half-integers to index irreducible representations. With this convention, the local plaquette ampli- 
tude ([2]) becomes: 

(14) Ap(/,p) = c,^(/3) = Up + l)e-wi''(='^+^\ 

By inspection of the plaquette amplitude, it is immediately clear that one can write down the plaquette 
dependent part of the interaction factor in closed form: 



n 
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where 5p{s) is the change in the spin labelling plaquette p due to the presence of the worldsheet s. The 
jp are of course a function of the background spin foam /; we suppress the / dependence for compactness. 
Because the interaction factor for a plaquette is simply one if it is not in the worldsheet s, we can rewrite 
the above as product of factors over the worldsheet (rather than the entire lattice) as follows: 

pes V ■jp + ^ J 

We discuss how the different factors in ([TB|) behave in both the self-avoiding and self-intersecting cases in 
Section 4. 



3.2. Vertex Interaction. We now turn to the part of the interaction factor not contained in the plaquette 
interaction Ip. Using a particular recoupling discussed in [5] and previously reported in [2Hl]j the edge 
amplitudes Ae and Ae can be set to unity and all the rest of the vacuum amplitude expressed as a product 
of 6j symbols of SU (2) associated to the vertices of the lattice. Appropriate modifications are to be made 
along the Wilson line F; we shall presently focus on the situation on the worldsheet but awaj0 from F. With 
this choice of recoupling, the interaction associated with the vertex can be written as 

5 
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In the above expression, the ai{v, /), . . . , fi{v, f) are the arguments to the product of five (indexed by i) 6j 
symbols associated with the vacuum spin foam / at vertex v. Upon introducing a worldsheet s, some of 
vertices will have their arguments displaced by some amount which we shall denote by Sl^{v, s), . . . , ^}(w, s). 
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^ The modified amplitudes along F will provide an interaetion factor that scales as the perimeter of F, whereas the rest of 
the worldsheet will be a product of a number of interaction factors that are bounded below by the area of F, and thus more 
relevant for confining effects. 
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3.2.1. Vertex interactions at large spin. We recall a well known approximation to the Racah-Wigner 6j 
symbol that has been conjectured to hold for large spins: 
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is the Cayley-Menger expression for the volume of the associated tetrahedron, which has side lengths given by 
the spin labels; the determinant that appears is the Cayley-Menger determinant. The formula assumes the 
limit is taken through arguments satisfying the Euclidean condition (the tetrahedron can be isometrically 
embedded in flat Euclidean space, a condition equivalent to the Cayley-Menger determinant has positive 
sign). 

Although suggested by Ponzano and Regge as early as [22], to the authors knowledge, the above asymptotic 
formula has been proven to hold rigorously only recently [T^[^51[^ and only in the specific case of a rescaling 
limit in which the arguments of some admissible 6j symbol are all scaled by the same integer k. In this case 
of rescaling, it was also shown that the non- Euclidean cases decay exponentially [12l[24]. Below we shall 
assume that the oscillatory or closely related exponential (a good discussion can be found in [M]) form of the 
Ponzano-Regge formula holds for any 6j symbol containing arguments of sufficiently large spin. This is in 
the spirit of the original conjecture but a slightly stronger result (stated in Section 5) than that of rescaling 
limits proven to to date. To see this, consider that at any spin scale there are tetrahedra with very long side 
lengths (high spin) which are not the uniform rescaling of any tetrahedra of smaller side lengths. 

Consider ly, which we shall define to be the contribution from a single vertex in the asymptotic limit: 
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Here V and 9'^ are the volume and set of dihedral angles associated with the tetrahedron whose side lengths 
have been displaced by the Si . We now consider the perturbed asymptotic formula for the 6j in more detail: 
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For concreteness and to motivate more general conjectures, consider the case of self-avoiding worldsheets, 
in which (5 = 1 (or (5 = 0, as not all spins local to vertex will be perturbed). In this case, changes in the 
dihedral angles b0i will go to zero in the limit of small perturbations about large "tetrahedra" . This leaves 
the term 'Y^^iib0i. It may not be immediately clear that this quantity is small since (although the b0i are 
getting smaller) we are studying the large ji limit. However, when ji are large, if E0i is close to zero, we can 
employ the fact that 

(22) \Y.^-^^^-^- 

i 
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due to the Schlafli identitji [TC]. We use this fact in the foUowing: for a self-avoiding worldsheet, we have 
Si{v, s) = 1 and thus find 



^^^^ V^cos [E.iMv) + 1)1 + (1 + S,{v, s)y-^+Mv)^ + f ] ^ ^ 

Similarly, since we are considering a ratio of V with the volume V' of the tetrahedron that has been perturbed 
by single units of spin, 

(24) ^ ^ 1 

in the limit of large ji . Combining volume and cosine factors in the asymptotic formula allows us to conclude 
that 

(25) /y(/,s)-l 

in the limit of large ji for self-avoiding worldsheets. We shall refer to this as the vertex decoupling property. 
While we have worked with the Euclidean case, the case of small perturbations in spin of the non-Euclidean 
follows by a similar analysis. 



4. Confinement in the weak coupling limit 

We now proceed to apply the general confinement scenario proposed in [7] to our present case, the first 
step of which separates contributions to the expectation value coming from asymptotic and sub-asymptotic 
backgrounds. 

4.1. Decomposition into asymptotic and sub-asymptotic backgrounds. Let a positive integer j* 
be given. It will be used to control the transition to asymptotic behavior for the vertex amplitudes and 
interaction factors. Background foams containing any plaquette spin less than j* are said to be sub-asymptotic 
and the set of such foams is denoted by j7o; the complement of this set is denoted by J'^ and its elements are 
described as asymptotic background foams. With these definitions, we can write the Wilson loop expectation 
value as follows: 

(26) (OrJ - (OrJo + (OrJ + 

where we have decomposed (Or„) into asymptotic and sub-asymptotic contributions, denoted by (Or„)+ 
and (Or„)o respectively. This allows us to apply different approaches to the two cases; we shall discuss each 
in turn. 

4.2. Asymptotic backgrounds. We decompose the first term of equation (I26p into the set of worldsheets 
that are self-avoiding, 5*^ and the remaining self- intersecting worldsheets, denoted by S'p. 

(27) (OrJ+ - Zy' E^"'^^(/'^)^^(/'^W) 

= E E + E E ^"'^^(/'^)^^(/'^)^(/)- 

We now consider the two types of worldsheets. 



similar argument was used by Roberts in |23l in relating asymptotic formulae for tetrahedra differing by one unit of 

length. 
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4.2.1. Self-avoiding worldsheets. We start with the case of self-avoiding woridsheets, as the situation is 
particularly clear in this case. Consider the interaction factor: 

(28) Ivif, s) = ^^evUUaAf) + S.As)}_ 

Ylvev 

Based on vertex decoupling noted in Section 3 equation ((25l) . we apply vertex decoupling to the ratio 
I^f = {^'■"(■^)+^'-"(^)} and find that it is bounded by 1 + efi*), where e ^ as 7* ^ oo. 
We next recall from Section 3 the form of the plaquette interaction amplitude Ip: 

As we are in the case of self-avoiding worldsheets, we have Sp{s) — 1 for all plaquettcs in s and can simplify 
to 

(30) /p(/,s) = []e-/3-'0/ + i) r^P±|') ^[]e-/5-^(w + l)(l + 0(£)) as j* ^ 00 

pes ^''^ ^ pes 

The notation 0(e) denotes a quantity that is no greater than e in magnitude, where e is a bouncQ that 
approaches zero as j* — > 00. Combining the vertex and plaquette interactions just discussed, let us now 
consider a general contribution to (Or^)^"", the self-avoiding part of (Or„)_|_: 

(31) (Orjf = V ^ ^ i?-i/p(/,s)/y(/,5)A(/) 

fej+ sGS^ 

where P(s) and V{s) denote the sets of plaquettes and vertices; absolute value lines around finite sets denotes 
the number of elements. To further resolve the factors that may contribute to or diminish confinement we 
organized the sum over worldsheets according to their area: 

(32) (Or J? = ^y' E E E i^-^e^"(«^(^)+5:,eP(., + o(£))5|v(.)I+|pMU(j) 

^ VE E E D~^e^~'^^^<^ns^^^)A{f) asj*,/?^^. 
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where ^Min is the area (in number of elementary plaquettes) of a minimal surface bounding F, and S^{A) 
is the restriction of the set Sp to surfaces of area A. 



We have written (|32p in a way that is particularly transparent with regard to confinement in both the 
strong and weak coupling limits. First, we see from the first line of (j32|) that in the strong coupling limit of 
/3 — » 0, area law behavior will eventually dominate any roughening effects — while the growth in the number 
of surfaces of a given area is exponential in the area, it is depends only on enumerative geometry and is thus 
clearly independent of /3; conversely, we are free to set (3 arbitrarily small and thus force a tension — this is 
has long been a standard result in the literature. 

By contrast, in the /3 — > cxd limit we find that the tension provided by the factor decays to zero (for 

any fixed 5;), and thus can not provide tension to general worldsheets. The determining factor thus becomes 
the sum over jp over the worldsheet surface, as can be seen from the second line of ((5^ . Roughly speaking, 
if the background spin averaged over any given worldsheet grows faster than /3 (which divides the sum) as 
/3 ^ 00, than a positive tension bounded away from zero will be present in the weak coupling limit. This 
particular means of realizing area damping via the plaquette interaction is presented as a formal conjecture 
in Section 5. 



As e appears in more than one type of bound each of which requires a certain minimum j * , we take the maximum of the 
j* for each type of bound to control the overall limit. 
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4.2.2. Self-interscicti'iig worldsheets. The set of r-admissiblc spin foams contains contributions whose reso- 
lution into background and worldsheet require a worldsheet that intersects itself. In general, intersection 
can happen along vertices, edges, and plaquettes (or any combination of these). We first consider a self- 
intersection that occurs at one or more vertices or edges, but not a plaquette. The important consideration 
here is that there are a limited number of ways this can happen without also changing a plaquette spin. 
Indeed, on a three-dimensional cubic lattice only one self-intersection can occur along an edge without in- 
creasing any plaquette spins; similarly for vertices. Because this sort of self-intersection will generically lead 
to 6p{s) = 1 in the vertex amplitudes of the intersection, we invoke the vertex decoupling shown above to 
find that these types of self-intersections do not provide any negative tension in the weak coupling limit. 

The remaining cases consist of worldsheets in which a self-intersection occurs at a plaquette. For some 
range of self-intersection {6i > 1), the same arguments for the self-avoiding case can be applied. Moreover, 
as j* gets larger this range will grow as it will take larger perturbations to appreciably change the angles 
in the asymptotic formula. However, at any j* there will be some sufficiently large Si such that the change 
in dihedral angles is too great to be within the assumptions of the self-avoiding case. What is required to 
deal with these cases is a condition that bounds the growth of Iv for increasing Sp{s). One way this could 
be true (described in more detail in Section 5) is if there existed a bound on ly that grew more slowly than 
the exponential decay factor that comes from Ip. If this is the case then the cut-off j* could be raised until 
the 6i appearing in the Type 1 factor from the Ip were large enough to exceed the bound on 7y, resulting 
in tension for these sheets. 

4.3. Sub-asymptotic backgrounds. Recall that for given a spin foam /, a sub-asymptotic plaquette is 
one for which jp{f) < j*, where j* is a fixed cut-off. We further define a spin foam to be sub-asymptotic 
if it contains one or more plaquettes that are sub-asymptotic. Our approach in treating this case will be to 
find a bound for (Or„)o that is suppressed exponentially in the weak coupling limit. 

We start by dividing the sub- asymptotic background spin foams into classes F'^{N), each element of which 
is a pair consisting of a subset of the lattice with A'^ plaquettes and an associated locally admissible spin 
foam labelling of those N plaquettes by sub-asymptotic spins and edges of these plaquettes by compatible 
intertwiners. We next split the sum into asymptotic and sub-asymptotic components, for both background 
spin foams and worldsheets, as follows: 

(33) (OrJo = Zy'Y. J2D~'lpif^^)W,s)A{f,s) 

/e Jo seSr 

= E E E E D-\f,r,s^,s*) 
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X Iv{f,8°)Iv{f\8*)Ip{f,s'')Ip{f\s*)A{f)A{r), 

where the set F*{p) consists of asymptotic completions of the sub-asymptotic foam these are assignments 
of spins to the set of asymptotic plaquettes and edge intertwiners compatible with a given assignment of spins 
to sub-asymptotic plaquettes and edges, /°. Given an /°, we define another set S^if^) that is the restriction 
of the set St to the sub-asymptotic plaquettes of Finally, for any given s° and the set of admissible 
completions of the worldsheet in asymptotic regions can be formed; we denote it by 5p(s'',/°). Note that 
lower dimensional simplices at boundaries (edges and vertices with both asymptotic and sub-asymptotic 
incidcmt plaquettes) between asymptotic and sub-asymptotic regions are defined to be in the sub-asymptotic 
component. 

A considerable obstacle to forming a useful boimd on (Or„)o is the fact that A{f) can take on nega- 
tive amplitudes. In particular, in forming upper bounds we will be careful not to take absolute values of 
amplitudes being summed over asymptotic regions of the state space, as quantities of this type can grow 
exponentially in the volume of the lattice faster than Zy itself (this manifests itself in numerical simulation 
as the decay of the sign expectation value), where cancellation occurs between positive and negative ampli- 
tudes. In contrast, for sub-asymptotic quantities at a fixed cut-off, we shall see that there arises factors that 
are bounded by exponential factors in the volume at a rate that is /? independent. We will show how they 
can be suppressed by terms coming from asymptotic factors in the (3 ^ oo limit. In summary, in proceeding 
we will take absolute values on sub-asymptotic quantities but not with asymptotic quantities. 
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We will start by considering the sub-asymptotic interactions and amplitude that appears in ([55]) : Iv{f^, s°), 
/p(/°,s°), and A(/°,s''). Recall that /y(/°,s°) can be written as a product of five factors of the form 

'^■^"l^-'^'^'^l^-'^ , a ratio of a 6j symbol^ with background spins displaced by the worldsheet and the original 
un-displaced 6j symbol. In the case where the 6j symbol contains all arguments less than j* , than the 
interaction associated with the ratio fa^'tj^^+'^t^)^ is bounded by 

(34) A^i = |(min[{J}l )-'\, 

the inverse of the smallest value taken by the 6j symbol on the finite set of values satisfying j < j* for all 
arguments j (because the numerator is a 6j symbol it is simply bounded by unity — see Appendix A). 

In the case where 6j symbols contain a mixture of asymptotic and sub-asymptotic arguments, we assume 
(currently based on numerical experiments rather than rigorous proof) that there is a finite range of values 
Jo of the asymptotic arguments, beyond which the 6j symbol becomes monotonically decreasing, and hence 

gives a ratio of ^''"^P^f^y^ ^^ < 1. In which case we have 

(35) X2 = |(min [{J}] )-'\. 

where the minimization is over all sub-asymptotic subsets of arguments, and for each sub- asymptotic 
subset, a minimum over the bounded non-monotonic range of the remaining arguments Jq. Defining 
M = max[A1i, A^2], we thus have 

(36) |/y(/°,s")| <A^5.o(.°)^ 

where the function vo{s^) is the number of vertices with some incident plaquettes labelled by sub-asymptotic 
spins; M. is raised to the fifth power as there are that many 6j symbols in the vertex amplitude. We shall 
also make use of the following inequalities for the remaining sub-asymptotic factors: 

(37) |A(/0,s")| <Pi, \Ip{f,s")\<P2, 

which follows from the fact that 6j symbols are bounded by unity; as well, the plaquette amplitude by 
inspection has a maxima and is bounded from below in the region where it is not monotonically decreasing 
with displacements by Si. The maxima always exists because the exponential parts of Ip will eventually 
damp out the linear increasing dimension factor for any non-zero /3. Note that the Pi and P2 are and 
dependent but we suppress this for compactness in notation. 

We turn now to the asymptotic summations and factors. For a given choice of and s^, we consider the 
factor contributed by summing over all completions of the asymptotic part: 

(38) r(/",sO)^ E D~\f,r,s^,s*)iv{r,s*)ipir,s*)Ain. 

/•6-F*(/°)^*eSf(/o,sO) 

We shall refer to T(/*',s°) as the asymptotic tail. Using this definition and introducing the sub-asymptotic 
bounds above, we can substitute into ([33]) and thus bound (Or„)o ^-s follows: 

A'=|P| 

(39) (Or Jo < ^y' E E E A^'"''^^°H^'i^2)^"^^)r(/0, s"). 

We next exhibit T{f'^, s") in terms of an average over a restricted partition function, which will suggests a 
conjecture that is critical for sub-asymptotic suppression in the weak coupling limit. We define 

(40) Z*{f)^ Yl ^(/*)' 



The statements made of 6j symbols in this section (essentially boundedness on finite sets of arguments) are meant to 
include the modified 6j symbols at T as well as the original unmodified Qj. 
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which can be thought of as a restricted partition function that sums over aU foams which have a fixed 
sub- asymptotic component f'^. We now write 

UU r,fO 0, _ ... .0. T.f>eFHfO) D-Hf, r, s*)Iv{r, s*)lAr, s*)Ain 
(4ij ,s ) - Z, (J ) =^ — — 

CSO 

A=q(/o,sO) 

where a is the smallest area of the restrictions of a F-filling surfaces to the asymptotic region; in general it 
depends on the position of T relative to the asymptotic plaquettes (e.g. ^Min if T is contained entirely in a 
region of asymptotic plaquettes) and to how the sub-asymptotic worldsheet s'^ ends on asymptotic regions. 
The sum up to "cx)" denotes that self-intersecting worldsheets (of unbounded amount of total spin) are also 
included. In writing (|41|) we have assumed a generalized form of asymptotic confinement: summing over all 
completions on (a possibly multiply connected set of asymptotic plaquettes) can be organized by area class, 
each with its own (positive) tension that is bounded away from zero by some (t*). Given this assumption, 
the sum converges, as it is a sum of terms bounded by exponentials with negative arguments proportional 
to area A. 

Inserting ([1T|) into ([55)1 we have 

(42) (OrJo < "i:' E E M'^"»^^"\P,P,r^^^ 0-^^) £ Oie-^^). 

N=l fOeFO{N)s»eS°{fO) ^ ^ A=a(P,sO) 

Let any N and Wilson loop V be given. For foams of a given fixed sub-asymptotic component /° £ F° {N) , 
let us characterize the possible terms in the in the inner sum: 

(43) M'^°^'''HPiP2r^'^ fZ2^\ g 0{e-^^). 

We recall that the minimal total area (both asymptotic and sub-asymptotic) of a worldsheet must be at 
least the area of the minimal surface spanning T. We see by inspection that the sub- asymptotic worldsheet 
interaction factor is bounded by the number of vertices (and hence plaquettes) in the asymptotic region — 
that is, there is a factor of negative tension 

(44) M^''"'-''°\PiP2Y"'-''^ = ^po{s}log{PiP2)^log(M)5vo{s") ^ ^T{r)Ao(s°) 

where ^0(5") is the number of plaquettes or vertices (whichever is greater) occupied by the sub- asymptotic 
component of the worldsheet s° and T{j*) = log(PiP2) -I- 51og(A4). This bounding factor certainly does 
not provide confinement as it has precisely the opposite effect: the weight of worldsheets s'^ increases 
exponentially with its area. However, the factor — ^ ' plays a crucial countervailing role. We consider 
the conjecture that the growth rate of the restricted partition function is slower by a factor proportional to 
the number of sub-asymptotic plaquettes (these are essentially degrees of freedom that are frozen) Ao(/o), 
that is: 

Zv ~ 

based on the assumption that both restricted and unrestricted partition functions increase exponentially in 
(3 and in plaquette volume. Because Ao{s^) < Ao{f°), we can form a bound that is directly in terms of s° 
and we do so on the rightmost side of ([^5]) . We then apply this inequality to find: 

N=\P\ 00 

(46) (Or Jo < E E e-^'^-^o^/") ^ JUlAois^) ^ 0{e--^) 

N=l f«eF°{N) sO£S°ifO) A=a{f°,s'>) 

N=\P\ I 00 



Consider next the quantity parenthesized above, which we shall define as S{f^): 



oo 



(47) S{f) = J2 e(^(^'*)-^'^)^«(^") Oie-^"") 

sOeSO(/0) A=a{p,s°) 

OO oo 



with the last inequahty holding for sufficiently large (3\ the inequality is saturated when (3 is tuned so as to 
provide tension r equivalent to the asymptotic tension and becomes a strict upper bound for larger [3. Note 
this last bound is independent of /°. 

We shall now introduce a number of bounds to control the growth of the remaining summation over N 
and F°{N). First, the number of partitions of plaquettes into N sub- asymptotic and \P\ — N asymptotic 
groupings is simply ('^')- We also observe that for a given number TV of sub-asymptotic plaquettes, there 
are no more than (j*)^ labellings of the sub-asymptotic plaquettes (of course not all will be components of 
admissible spin foams). 

We now combine the various bounds introduced above to find an overall bound on (Or„)o- Recalling that 
N = Ao{f'^) is the number of sub-asymptotic plaquettes, we have 

Ar=|P| oo 

(48) (OrJo <P^oo E E e--'''"''''' E 

N=l /0gi?0(7v) A=AMi„ 

^ E ('^')(r)^e-^'^^ Y ^(-"^■') 

N=l ^ ^ A=AMi„ 
\P\ oo 

< Y e('°s(l^'l)+iog(r))JVg-f,3Ar Oie'"'^) 
where we've used the inequality (1^1) < \Pf = (,iog{\p\)N ^ -y^g 

assume that in the limit of large (3 invoked 
here, the lattice volumes (as measured by number of plaquettes \P\) that fall within the scaling window 
grow no faster than exponential in /3, and hence the theory has the appropriate critical behavior to provide 
a continuum limit. 

Finally, we state our conclusion for the sub-asymptotic case: subject to the hypothesis given in (f45|) for 
sufficiently large (3, sub-asymptotic contributions 

(1) Are hounded by a decaying exponential series with leading order having tension associated with 
asymptotic tension. This series converges. 

(2) The coefficient in front of the convergent series over worldsheet areas is a decaying exponential for 
sufficiently large (3, for any lattice volume (as measured by number of plaquettes). The rate of decay 
is also proportional to the number of sub-asymptotic plaquettes N . 

On account of point (2), we conclude that in the weak coupling limit, all suh- asymptotic contributions are 
suppressed: (Or,„)o ^-t a fixed cutoff j* (e.g. as given by some bound e) can be made arbitrarily small for 
sufficiently large (3. 

5. Conjectures 

In Section 4 we described scenarios in which area law behavior for asymptotic backgrounds could emerge, 
and how sub- asymptotic backgrounds are suppressed — both in the weak coupling limit. In this section we 
collect together more precise statements of those properties that are sufficient to provide area law behavior, 
as well as an additional balance condition, which if true would overcome the technical problems raised by 
the non-positivity of the spin foam amplitudes. 

Below we have presented these conjectures according to the general structure proposed in [7]. That is, we 
start by conjecturing that the growth rate of background spins is of the same or high order as that of (3 itself 



in the /3 ^ oo limit. We then introduce sub- asymptotic suppression to argue that asymptotic background 
foams are the only ones relevant in the weak coupling limit. Next we take up the question of vertex stability 
on asymptotic backgrounds, which breaks up into two cases — self-avoiding and self-intersecting worldsheets. 
In the self-avoiding case, we conjecture that the Ponzano-Regge asymptotics may be true in a uniform sense, 
which would (by the analysis of Section 3) be sufficient to establish vertex stability in the self-avoiding section 
(for all asymptotic spin foams). To treat the remaining self-intersecting sheets, we conjecture the existence 
of a lower bound on the 6j symbol (outside regions of monotonic decrease in Sp) that is sub-exponential in 
Sp (introduced in Section 3). Such a bound would be sufficient to provide vertex stability, as the relevant 
part of the plaquette interaction decays exponentially in dp. 

Given the preceding conjectures, our final conjecture is that the balance condition introduced in [7] is 
true, allowing us to form an area-law upper found on the growth the Wilson loop expectation value. 

5.1. Area damping from plaquette interaction. Considering both avoiding and self-avoiding sheets, we 
define the plaquette tension as follows: 

(49) /p(/,s) = (^c'^'' SJs):}^"^ (^c^t^^)'' T.^..SJs)(SJs)+2)'j I -Q jp + 1 + Spjs) j 

Vpes ~^ J 

and apply the definition of area damping from plaquette interaction introduced in [7] 



(50) 



EfeF. \EA>A^,^.iT)EseSriA)D-Hf,s)e-^^(f'^)^\ A(/) 



< oo as P 



Define the first, second, and third parenthesized factors of (|49|) to be Type 1, 2, and 3 factors, respectively. 
At any fixed (3 and background foam /, the Type 2 factor will eventually damp out growth in the Type 3 
factor. However, as, [3 ^ oo, the number of terms in the Type 3 factor that could give above unity overall 
Ip grows without bound due to the weakening of the Type 2 factor with inverse (3. 

By elimination, we conjecture that the Type 1 factor is critical to obtaining an area law. We observe 
that on backgrounds in which the average spins over any sheet is greater than (3, the first factor will provide 
a positive tension to all worldsheets, with self-intersecting ones being suppressed at a rate proportional to 
their 5{s). As the asymptotic cut-off scale is increased, the slope of the linear growth in 5 in the Type 3 
factor will insufficient to overcome the tension from Factor 1. Thus for sufficiently large cut-off it is possible 
for Ip to provide a sub-unity (positive tension factor) in the /? — > oo limit for a large class of background 
spin foams. Thus, as long as backgrounds of this type tend to dominate contributions to the expectation 
value in the /3 — > oo limit, the area damping condition will be met. More succinctly: 

Conjecture 5.1. As (3 oo, the set of background spin foams for which the Type 1 factor provides a 
positive tension sufficient to satisfy / f50|] dominates the contributions to (Or„)- 

This is simply a more specific form of the general area damping condition given in [7 , as it identifies the 
means by which the defining inequality ()50p is realized. We remark in passing that this conjecture is well- 
suited to numerical simulation — the expectation value of the Type 1 factor on the vacuum ensemble can 
be computed using both conventional and dual spin foam algorithms (as in [S]) and is the topic of ongoing 
investigation by the author. 

5.2. Sub-asymptotic suppression. Let Z*{f^) be defined as in equation (|40p . We conjecture that 
Conjecture 5.2. For any sub- asymptotic spin foam component f'^ we have 
(51) < e-/^Mo(/°) 

Zy 

where Ao(/'^) is the number of plaquettes in f'^ and pi is a positive constant. 

The number of degrees of freedom summed over in forming Zy and its restrictions Z*{f'^) is proportional 
to the volume of the complement of /°. 

Testing this conjecture is possible by dual spin foam methods [5| by studying the characteristic observables 
of /" — functions that are unity on background foams / where spins are consistent with /o and zero otherwise. 
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It would also be interesting to investigate analytic or numerical methods of showing the conjecture that used 
the integral presentation of the partition function as is done in conventional lattice gauge theory. 



5.3. Vertex stability: self-avoiding worldsheets. The following conjecture allows us to use an asymp- 
totic form of the vertex amplitude for large values of spin j. 

Conjecture 5.3 (Asymptotic formula with uniform cut-off). Let e > be given. Then there exists some j* 
such that if ji > j* for all of the ji ■ ■ -je arguments to the 6j symbol {jv} then 

(52) -%t- = l + e(r) 

Asy(j„) 

where 

(53) Asy(j' ) = I \/ 3^ '^0^ + 1)% + f } f°'^ Euclidean tetrahedra. 

[ Otherwise, decaying exponentially, 
where the Euclidean condition and other aspects of the asymptotic formula are discussed in Section 3. 

This allows one to analyze the asymptotic case of self-avoiding worldsheets as in Section 4.2.1, culminating 
in ([5^ . at which point vertex stability (for that class of worldsheet) reduces to the area damping conjecture. 

As mentioned above (Section [3.2. ip . rigorous derivations of the asymptotic formula in the literature are 
presently limited to rescaling limits, where a set of admissible arguments is rescaled. Recent advances in 
analyzing 6j symbols [nil2j and in particular [14j are encouraging — we expect it may be known soon rather 
the Ponzano-Regge formula or some generalization thereof can provide a uniform bound that behaves in the 
way stated. 

We should note here the significance of the vertex interaction going to unity (rather than above or 
below unity, which would provide negative and positive tension factors). One implication is that for self- 
avoiding worldsheets on asymptotic backgrounds, confinement must by elimination come from the plaquette 
interaction. Further, assuming confinement is true, this result can be used as a starting point (base case) in 
showing that the self-intersecting case also has decoupled or at least provides bounded negative tension. We 
turn to this case next. 

5.4. Vertex stability: self-intersecting case. In [7] a relatively "high-level" definition of vertex stability 
is given, with the idea that different types of estimates may be needed for models of different gauge groups 
and dimensionality. An area averaged tension is defined as 



(54) e 



.r(f,A)A ^ E^6Sr(A)_^v(/,s)/p(/,s) 



gTo(A-AMin) 

this allows comparison of the interaction factor with the "negative tension" arising from the exponential 
growth in surface number, n{A) — e'^f^'^"'^"'"), leading to the condition that 

\p L-To^Min V p{To-f{f,A))A] A( f\ 

(55) (Or„) = 



EfeF.Aif) 
indicating that 

(56) f{f,A)-To>0 

is the characteristic inequality. Showing that this inequality is satisfied such that the infinum of f (/, A) over 
all / and A that aren't suppressed as f3 ^ oo then defines the condition of vertex stability for a model. 

A general approach to establishing the above is to carry out a local analysis, in the sense that if |/v(/, s)\ 
is above unity over a worldsheet, then some factor of which it is a product must be over unity. Thus, if one 
considers all admissible arguments to a factor of |/\/(/, s)| at a single vertex and establishes bounds on that 
quantity for all admissible arguments, one could extend it to a bound over any worldsheet for any asymptotic 
background. 

In our present case we can pursue this local strategy, and moreover exploit the fact that vertex interaction 
factors are directly related to perturbations of 6j symbols. The availability of asymptotic estimates then 
provides a means of realizing the necessary condition, leading us to conjecture 



Conjecture 5.4. We define sub-exponential in Si to mean that exj^T^^'^i )) ~* ^ '^^V ^' '^'^^ 
be any collection of the Si arguments appearing in |i7p . For all possible admissible arguments where ly is 
non-zero, one of the two is true: 

(1) Monotonic decreasing with Si or 

(2) \Iv{fTs)\ is sub- exponential in Si. 

The expectation behind this conjecture is that in the Euchdean regime where the 6j asymptotic is charac- 
teristically oscillatory (multiplied an effective power law volume factor) there will exist some sub-exponential 
bound away from zero. This will in turn limit the rate at which displacements of Si can lead to growth rela- 
tive to undisplaced 6j symbol; if so, then the exponential damping from the Type 1 plaquette factor will be 
able to attenuate worldsheets with arbitrarily high levels of self-intersection. 

In the non-Euclidean cases where the Qj may exhibit exponential decay in the sum over jiOi, large dis- 
placements of the 6j factors in |Jy| will be less than unity if 6j is monotonic decreasing. Monotonic decay 
when increasing single arguments or some subset of the arguments has (to the author's knowledge) not yet 
been rigorously shown, but is well-supported by numerical testing by the author and in the literature. 

5.5. Balance relations among positive and negative sign contributions. As emphasized in 7 , the 
amplitude A{f) is in general non-negative for dual spin foam models of Yang-Mills theories and that is 
certainly true of the present case; this raises some considerable technical complications in forming an upper 
bound. We next make the following conjecture: 

Conjecture 5.5. The interaction factor for any given worldsheet s, averaged over negative amplitude states 
Fy approaches its average over positive amplitude states Fy as /3 — > oo. 

What this allows us to do is take either the positive or negative weight background spin foams to estimate 
the Wilson loop expectation value. Taking the positive weight configurations: 

V .p-T^oAMinY- (n-f(f,A))A\A(f\\ 

^''^ ^""'^^ ^ E,... 

< (--"°"-i:.,...„e--^)E,..,..HA(/)i 
" E,,.r-i^(/)i 



e 

A>A 



where the bounding tension is defined as r* — inf^T*+ ^ [t(/, A) — Tq]. The term £'(/?) represents contribu- 
tions where the existence of a bounding tension does not hold, but which can be shown to go to zero in the 
f3 —>■ oo limit. 

An alternative to the balance condition is to conjecture that there is a subset of J-gp of stationary phase 
which dominates the partition function in the continuum limit. In this case, assuming numerator of the 
expectation value (as a sum of A{f) weighted by the worldsheet observable) is also dominated by config- 
urations from J-^sp, then we can use positivity on the set J-gp to form a similar bound. The availability of 
asymptotic estimates for Ay combined with the suppression of sub-asymptotic configurations may render a 
stationary phase tractable, and is the subject of ongoing work by the author. 



5.6. Final steps. The conjecture (|57)) provides a presentation of the expectation value for a single fixed 
Wilson loop r as a decaying exponential sum over classes of worldsheets of increasing area. To consider 
confinement, one must consider the behavior across a range of Wilson loops corresponding to increasing 
spatial separation (and thus increased area) of quark- antiquark pairs. As discussed in [7], one way to 
provide this is to define 

(58) fe=infg[T*(r)]. 

The (leading order) of each world-sheet area expansion of (Or„ } can thus be bounded by a function propor- 
tional to e~'^ -^Min^ |;]^yg yielding the area-law decay characteristic of a confining theory. It should be noted 
that such a bound does not (by itself) rule out local variations in the expectation value that may fluctuate 
within the exponentially decaying envelope. 
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We mention in passing how a similar analysis could be applied to the mass gap observable (2-point 
correlation function), and in principle higher order correlation functions. As discussed in Appendix B of [7], 
the main point of departure is that the topology of worldsheets is different — the background foams are 
unchanged however the worldsheets terminate on two fundamental plaquettes rather than single loop F. As 
a result, the rate of growth with increasing area may expected to differ somewhat, as would the effective 
tension of a given area class. 

6. Conclusions 

As described in [7] , a Wilson loop expectation value in conventional lattice gauge theory can be expressed 
as the expectation value of a certain type of observable on the vacuum spin foam state sum. This observable 
is a sum of interaction factors — one for each worldsheet bounding F (modulo a degeneracy factor). Each 
interaction factor is manifestly a product of factors over plaquettes and vertices of the worldsheet. Due 
to the factorization of amplitudes over worldsheets and the ability to organize worldsheets by increasing 
area, this spin foam formalism is particularly well-suited for identifying the absence or presence of area law 
decay. In the present paper, we introduced the interaction factors associated to the choices G = SU{2) and 
space-time dimension Z) = 3 specifically and described methods for understanding their behavior in different 
sectors: sub-asymptotic, asymptotic with self-intersecting worldsheets, and asymptotic with self-avoiding 
worldsheets. 

The resulting analysis shares a number of similarities with the study of lattice gauge theory through strong 
coupling diagrams in the strong-coupling limit. The weak coupling case considered here however requires a 
priori consideration of all vacuum spin foam backgrounds as well as all worldsheets bounding F, including 
those with arbitrary amounts of self-intersection. By contrast, the theory at strong coupling depends on a 
much more restricted class of diagrams. 

The scenario analyzed here follows the general ansatz described in [7], where background foams that 
dominate the weak coupling limit provide tension to worldsheets through plaquette interaction factors. 
In order for this scenario to hold true, the stability of vertex interaction property needs to be shown. A 
remarkable result of the present work is that vertex stability for self- avoiding worldsheets follows immediately 
from a well-known asymptotic formula for 6j symbols, of which the vertex interaction is a ratio of products. 
For worldsheets of arbitrary self-intersection, we explain how the vertex interaction factor (and thus, the 6j 
symbols that comprise it) can be bounded in such a way that the plaquette interaction can still provide a 
tension. 

We further show that subject to a conjecture on the relative growth rate of restricted sums relative to the 
full partition function, backgrounds with sub- asymptotic spins (below a given cut-off) can be ignored in the 
weak coupling limit. This allows the remaining conjectures to be described in terms of the asymptotic form 
of the vertex amplitudes, which should be susceptible to a variety of analytic methods. 

Despite the above results, the picture currently presented by this approach is incomplete in several re- 
spects. The behavior of local amplitudes critical for vertex stability involves several properties relating to 
the asymptotic behavior of 6j symbols which, while suggested anecdotally by numerical experiment, remains 
to be proven rigorously for all admissible values. At the level of the statistics of the model, the growth 
rate of background spins relative to (3 remains to be shown. Further statistical results are also require to 
assess the growth rate of restricted relative to unrestricted partition functions, the essential ingredient of 
sub-asymptotic suppression. Finally, the presence of alternating signs continues to cloud the picture, and it 
will be interesting to see if stationary phase methods may shed light on this difhculty. 

In summary, the present work gives a working example of how spin foams can be used to arrive at a 
gauge-invariant understanding of confinement at weak coupling similar in concept to the standard proof 
of confinement at strong coupling. To provide a concrete scenario for confinement, a number of specific 
properties of the spin foam model with varying levels of heuristic and numerical support were conjectured 
which taken together are sufficient for an area law bound. We hasten to add that there could be a number 
of alternative realizations should one of these conjectures fail — nonetheless, we expect that the general 
framework we have illustrated will increasingly constrain the behavior local amplitudes and statistical be- 
havior of the model as more results (both negative and positive) are obtained, to the point where alternative 
possibilities can be systematically explored. We do hope the results obtained thus far are sufficient to spark 
further progress on the conjectures given here as well as motivate similar investigations in four dimensions 
and for G = SU{3). 
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A. Absolute Convergence of vacuum and F-admissible spin foam state-sum 

Theorem A.l. (Absolute convergence of Z). The partition function Z||, in which every the absolute 
value of every summand in the original partition function is summed, exists for any (finite) lattice k. 

The Racah-Wigner 6j symbols used in the present work are bounded by unity, shown in P^llll] for example. 
Therefore we can write: 



(59) 



feF^ \i€Uf 



In the final step we replace a sum of products by a product of sums. Each sum in the product can in turn be 
bounded by an integral, whose existence is guaranteed by the finiteness of integrating a polynomial against 
a gaussian on M+ : 




(60) Z||(/3,«;) < 1 M 2J Uf + l)e-m^f^^ /+^^^' < / djfiff + l)e-^^f^Jf+'' < oo 



where \F\ are the number of plaquettes (polygonal faces) in the 2-complex k. 

The proof goes through essentially without change in the case of the charged partition function, except in 
this case both 6j symbols and modified 6j symbols that arise in vertex amplitudes along T (see 'W) have to 
be shown to be bounded from above. The modified 6j is a special case of a unitary evaluation of an SU{2) 
spin network (without loops or trivial components), and thus a constant upper bound is provided by the 
work of Abdesselam |1, , which follows on a conjecture given in fTT by Garoufalidis and Van der Veen. 
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